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We propose an ultra- thin elastic cloak to control the scattering of bending waves in isotropic heterogeneous 
thin plates. The cloak design makes use of the scattering cancellation technique applied, for the first time, to 
the biharmonic operator describing the propagation of bending waves in thin plates. We first analyze 
scattering from hard and soft cylindrical objects in the quasistatic limit, then we prove that the scattering of 
bending waves from an object in the near and far-field regions can be suppressed significantly by covering it 
with a suitably designed coating. Beyond camouflaging, these findings may have potential applications in 
protection of buildings from earthquakes and isolating structures from vibrations in the motor vehicle 
industry. 

I n the past few years, metamaterials have become an attractive research focus for physicists and engineers due to 
I their ability to control electromagnetic, acoustic, and elastic waves in ways that cannot be achieved using 
I "regular" materials. Consequently, exciting applications of metamaterials ranging from sub-wavelength 
imaging to cloaking have been reported 1 " 22 . Soon after the advent of negative refraction and sub -wavelength 
imaging 1 , it has been shown that dielectric or conducting objects could be made "invisible" to electromagnetic 
waves by coating them with plasmonic materials 2 ' 3 . A promising route to electromagnetic invisibility relies on the 
scattering cancellation phenomenon, where the fields scattered from a coating with negative (local) polarizability 
are specifically designed to cancel the fields scattered from the object 2 ' 4 " 7 . The scattering cancellation technique 
has been proven relatively robust to changes in the design parameters, geometry, losses, and frequency of 
operation 8 9 . It has also been used to make electrodes or metal films transparent to THz 10 or GHz 11 radiation. 
Experimental evidence of such invisibility devices has recently been demonstrated in the microwave regime 12 " 14 . 
An alternative route to designing invisibility devices makes use of coordinate transformations 1516 to mold the 
wave flow in desired manners 17 " 22 . It should also be mentioned here that with their possible applications in stealth 
technology, noninvasive probing, and sensing, electromagnetic invisibility cloaks open promising avenues for 
new technologies in medicine, defense, and telecommunications 23 ' 24 . 

Cloaking ideas and designs for acoustic waves have also been developed. The possibility of two-dimensional 
acoustic cloaking for pressure waves in transversely anisotropic fluids has recently been shown 25 . Subsequent 
investigation of this cloaking idea for concentric layers of solid lattices behaving as artificially anisotropic fluids in 
the homogenization limit has been provided 26 . Following a similar approach, a micro -structured metallic struc- 
ture has been designed as a cloak for surface liquid waves; the design has been validated experimentally around 
10 Hz 27 . Additionally, three-dimensional acoustic cloaking for pressure waves in fluids has been envisaged; this 
cloak makes use of the fact that the scalar wave equation retains its form under geometric transforms 28 ' 29 . Like in 
electromagnetic cloaking, the scattering cancellation technique offers an alternative path towards acoustic cloak- 
ing 30-32 . However, cloaking fully coupled pressure and shear elastodynamic waves in solids is a more elusive task, 
since the Navier-Stoke's equations do not retain their form under geometric transforms 33 ' 34 . Although the 
theoretical foundation of these cloaking effects is well established 35 ' 36 , the inherent requirement of singular bulk 
metamaterial properties that can hardly be met in practice still hinders the practicality of cloaking devices 37 ' 38 . 

In this Report, a biharmonic wave equation with appropriate boundary conditions, which describes the 
propagation of bending waves, is derived from the generalized elasticity theory 39,40 and is used in designing a 
cloak that operates in thin plates. The idea behind the cloak design is to choose a coating material with density p c 
satisfying p c < p 0 or p c < 0 and to conceal an object with density p s > p 0 . Here, p 0 represents the density of the 
surrounding medium. 
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The parameters of the cloak are deduced using the scattering 
cancellation technique applied to a platonic metamaterial coating 2 . 
Note that the term "platonic" was coined by McPhedran in reference 
to thin plates while investigating wave propagation in complex thin 
plate elastic artificial materials 41 " 43 . Recent advances in the field of 
high frequency homogenization in platonics provide a theoretical 
foundation for understanding frequency dependent effective mass 
density taking negative or less than unity values 44 . 

We then analyze the response of the platonic cloak concealing a 
cylindrical obstacle in the presence of a plane wave elastic excitation 
(harmonic vibration of the plate in the vertical z-direction). It is 
assumed that the out-of-plane dimension of the obstacle is negligible 
compared to its in-plane dimensions 41,42 . We show that in the quasi- 
static limit, i.e. for /c 0 a 5 ^Cl, where k 0 is the bending wavenumber in 
the surrounding medium and a s is the in -plane dimension of the 
scatterer, the scattering is dominated by the zeroth-order multipole, 
unlike in the electrodynamics case where the first significant order is 
the dipolar one. This is not the only fundamental difference between 
electrodynamic/acoustic and elastic wave scattering phenomena: the 
fourth order biharmonic partial differential equation, which typically 
describes the propagation of bending waves in ultra-thin plates, is not 
equivalent to the vector/scalar wave equations that describe electro- 
magnetic or acoustic wave propagation. For instance note that prop- 
agating and evanescent waves coexist even in homogeneous thin 
plates. Consequently, one can anticipate that a new cloaking method 
and new relevant physics are introduced following this route. 

Results 

It has been recently shown that the elasticity equations are not invari- 
ant under coordinate transformations 33 . Indeed, equations describ- 
ing the propagation of elastodynamic waves under a time harmonic 
dependence e~ imt read (in a weak sense) 

V-C : Vu + p 0 co 2 u = 0. (1) 

Here, p 0 is the (scalar) density of an isotropic heterogeneous elastic 
medium, C is the rank- four symmetric elasticity tensor (which has up 
to 3 4 entries even for thin plates), w is the wave angular frequency, 
and u is the associated three- component vector displacement field. In 
most practical applications, the lateral dimension is significantly 
larger than the plate thickness making the Kirchhoff assumption 
for thin plates applicable 39 ' 40 . Under this assumption, shear deforma- 
tion and rotary inertia are negligible. For such a structure, Eq. (1) 
reduces to the scalar fourth order biharmonic equation 

v\E(rfl 2 v[[p(r)/p,]- ll2 V-(E(rfl 2 VU(r,e))] } - 

b(r)/p o ] 1/2 **tf(r,0) = O, 

where U(r, 0) is the displacement in the vertical z-direction (assum- 
ing a cylindrical coordinate system due to the symmetry of the struc- 
ture), E(r) is the relative Young modulus, p(r) is the density, 
kl = co 2 p 0 h/Do, D 0 , and p 0 are the bending wavenumber (to the 
power four), rigidity and density in the thin plate, respectively, and 
h is the plate thickness. Consider the structure shown in the inset of 
the Fig. 1. An object of radius a s is coated with a shell of outer radius 
a c . For r < a s (inside the object), p(r) = p s and the bending wave- 
number is k s . For a s < r < a c (inside the shell), p(r) = p c and the 
bending wave number is k c . Object, shell, and thin plate all have the 
same rigidity D(r) = D 0 , same relative Young modulus E(r) = 1, and 
same Poisson's ratio v(r) = v 0 . Without loss of generality, we assume 
that this structure is illuminated by a plane wave propagating in the 
x- direction. The displacement field of the incident plane wave is 
expressed as U mc (r,6) = e lk ° rcose , equivalently it can be expanded as 

GO 

U mc (r,0)= s n i n J n (kor) cos n6, where the coefficients s 0 = I and 

n = 0 

s n = 2, n > 2. The scattered field, (7 sca (r, 6), must be finite at r = 0 and 
it satisfies the radiation condition at r — > oo. Thus, it is expressed as 




Figure 1 | Geometry under investigation: An object at the center of a thin 
elastic plate lies on the trajectory of a planar bending wave, incident from 
right to left in presence of a cylindrical shell that may act as a platonic 
cloak. Inset shows a top view of the cloak and object. 



GO 

U sc \r,0)= ^£ n i n [A n H%Xk 0 r) + B n K n (k 0 r)]cosn6, r>a c . (3) 

n = 0 

Similarly, fields inside the shell and the object are expressed as 

GO 

U**(r,ff) = K Y„(k c r) + D„K n (k c r) + EJ n (k c r) 

+ F n I n (k c r)\ cos n6, a s <r<a c , (4) 

GO 

U int (r,0)= J2 s ni"[GJ„(k s r) + H„I„(k s r)]cosn9, r<a s , 

n = 0 

respectively. Here, H^\.), /„(.)and !„(.), and 7„(.)and K n (.)are cylin- 
drical Hankel functions of the first kind, Bessel and modified Bessel 
functions, and Bessel and modified Bessel functions of the second 
kind, respectively. To solve for the coefficients in the above equa- 
tions, continuity of the field U, its derivative in radial direction d U/dr, 
the bending momentum M r , and the radial component of the gen- 
eralized Kirchhoff stress V r is enforced at the boundaries dXr = a s and 
r = a c , for each azimuthal order n (see supplemental materials 45 for 
explicit expressions of the bending momentum and the generalized 
Kirchhoff stress in cylindrical coordinates). This yields a matrix sys- 
tem of equations in scattering unknown coefficients A n and B n . The 
far-field scattering amplitude (or differential scattering cross-sec- 
tion) f(0) = V2?e- i{kor - n/4) lim U sc \r,0) is a measure of the 

r— >co 

cloaked object's visibility in direction 0 46 . The total scattering 
cross-section, cr sca , is the integral of f(6) over all angles, i.e. 

p27T 

(7 sca = 1/2 dO\f(0) | 2 . It may thus be expressed as 
Jo 

4 00 
*o»=o 

Note that coefficients B n are absent in Eq. (5) since the modified 
Bessel functions K n (k 0 r) have no contribution to the scattered field 
as r — > oo. Generally speaking, the possibility for an observer to detect 
the object in the far-field is determined by the value of cr sca . As a result, 
minimizing or completely canceling cr sca would lead to the undetect- 
ability (invisibility) of the object in the far-field, independent of the 
observer's position. This can be achieved by canceling the coefficients 
A n =A n / d n that significantly contribute to the scattering. Here, A n is 
given by the determinant: 
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-Jn(k 0 a c ) K n (k 0 a c ) -Y n (k c a c ) -K n (k c a c ) -J n (k c a c ) -I n (k c a c ) 0 0 

0 0 Y n {k c a s ) K n (k c a s ) ] n (k c a s ) I n (k c a s ) -J n (k s a s ) -I n (k s a s ) 

k 0 J' n (k 0 a c ) k 0 K' n (k 0 a c ) -k c Y' n (k c a c ) -k c K' n (k c a c ) -kj' n (k c a c ) -kj' n (k c a c ) 0 0 

0 0 k c Y' n (k c a s ) k c K' n (k c a s ) k c J' n (k c a s ) k c I' n (k c a s ) -kj' n (k s a s ) -kl' n (k s a s ) 

S Jn (k 0 a c ) S Kn (k 0 a c ) -Sy„(fc c a c ) -S Kn (k c a c ) S Jn (k c a c ) S ln (k c a c ) 0 0 

0 0 S Y „(k c a s ) S K „(k c a s ) S,Jk c a s ) S h (k c a s ) -S Jn (k s a s ) S h (k s a s ) 

T In (k 0 a c ) T, n (k 0 a c ) -T Yn (k c a c ) -T Kn (k c a c ) T In (k c a c ) T In (k c a c ) 0 0 

0 0 T Y (k c a s ) T K (k c a s ) T, (k c a s ) Ti (k c a s ) -Tj (k s a s ) Tj (k s a s ) 

(6) 

where 

S Zn (kr) = D 0 [n\l - v 0 )(/c z r) 2 ]Z n (/c z r) - D 0 (l - v 0 ) krZ^kr), 

(7) 

T Zn (k ir ) = D 0 [n\l - Vo)}Z n (kir) -D 0 [n\l- v 0 )(^r) 2 ] fc,r ^(fcr), 

where /c,e{/c s ,/c c? M and Z n {)e{J n {)J n {)M { n\M)J n {)}. 
Upper and lower signs in Eq. (7) should be selected for 
Z n (.)e{J n (.),Y n (.),H^(.)} an&Z n (.)e{K n Q,I n (.)}> respectively. The 
expression of d n could be obtained from Eq. (6) by replacing J n {.) by 
H n (.) in the first column (see supplemental materials 45 ). Given the 
general complexity of this expression, it is instructive to analyze the 
low frequency limit corresponding to elastically small obstacles and 
shells, i.e. k c a c <^l, k s a s <^.l, k 0 a c <^l, and k 0 a s <^l. Note that with the 
parameters in our study, it is sufficient to impose only /c 0 a 5 ^Cl and 
k 0 a c ^l. Under this assumption, cr sca in Eq. (5) is dominated by the 
monopole term A 0 : 

^~yJM 2 ~^^[p-P 0 +Ap- Pc )}\y=a s la c . (8) 

Note that here, terms scaling with (k 0 a c ) m , m > 3, in expressions of 
A n , n > 0, are safely assumed to be zero for k 0 a c < < 1. It is clear from 
Eq. (8) that cr sca of elastically small obstacles in a thin plate scales with 
k\a A c . Also, as y — > 1 and p c — > p 0 (no shell, only bare object), Eq. (8) 
reduces to 

By enforcing A 0 = 0, i.e. er sca ~ 0, in Eq. (8), we can derive the 
quasistatic design rule for cloaking: 



which relates p a p s , a s , and a c . It is found that, to satisfy Eq. (10), p c 
must take negative values for p s > p 0 and values larger than p 0 for 

Ps < Po- 

For larger objects, retardation effects become important and 
numerical calculations using the full 8X8 matrix in Eq. (6) are 
necessary to analyze the problem. 

It is interesting to investigate two additional types of obstacles, 
namely when a "clamped" and a stress-free boundary conditions are 
set on the inner boundary of the system (r = a s ), while the same 
conditions of continuity of (7, dUldr, M r and V r as before are hold on 
the outer boundary (r = a c ). For the first case, we have U = dU/dr = 0 
whereas M r and V r are unconstrained. The 8X8 system in Eq. (6) 
reduces thus to a 6 X 6 system (see supplementary materials 45 ). 
Applying the analysis described before to this case reveals that the 
monopole scattering coefficient A 0 ~ —1 in the quasistatic limit 
regardless of the cloak's parameters. This means that the resulting 
scattering amplitude/ (0) ~ — 2/ \fnk (or the scattering cross-section 
0-sca ^ 4/£) becomes singular in the quasistatic limit (all other scatter- 
ing coefficients tend to zero). This is unique to elastic waves, because 
a very small object would scatter infinitely and it has no equivalent in 
electromagnetics or acoustics 46 . This discussion shows that it is not 
possible to use shells of tailored density, in the quasistatic limit, to 
cancel the monopole scattering of a clamped obstacle since their 
corresponding coefficients differ by a factor (k 0 a c ) 2 . 



For the stress-free boundary conditions, M r and V r are set to zero 
on r = a s whereas the same continuity conditions hold for the other 
parameters. We obtain here also a 6 X 6 system (see supplementary 
materials 45 ). By applying the analysis above, we show that the mono- 
pole scattering coefficient of a stress-free hole of radius a s cloaked 
with a shell of density p c and radius a c = ajy, is given by 
— in 

A 0 ^ [l +y 2 (v + l)/(v - 1) + (y 2 - l)p c /p 0 ] . We see then that 

8 

it is possible to find values of the density of the shell p c and ratio y that 
make A 0 vanish thereby canceling the scattering from the cloaked 
object in the quasistatic regime. This condition is found to be 
PcM>=[l+y 2 (v+l)/(v-l)]/(l-y 2 ). 

In all the results presented in the rest of the Report, in-plane 
dimension of the object is set to a s = 1 m. Figure 2 shows the 
dependence of the scattering amplitude \f(0)\ on k 0 a s for p s /p 0 = 
0.1 [Figs. 2(a), 2(b) and 2(c)] and p s /p 0 = 10 [Figs. 2(d), 2(e) and 
2(f)]. We show here that curves of scattering amplitude \f(6) \ vs. 6 in 
logarithmic scale contain angular information about the scatterer 
[Figs. 2(a), 2(b) and 2(d), 2(e)] and curves of normalized backward 
scattering amplitude [f(7r)|/v^s vs - ko a s [Figs. 2(c) and 2(f)] give 
insight into the spectral dependence of the object visibility for an 
observer placed at the angular position of the illumination (mono- 
static scattering). The simulations are carried out for the first scen- 
ario (soft object e.g. soil) for k 0 a s = 0.1 (quasistatic limit) and k 0 a s = 
1 (Mie scattering). It is clearly seen that, for the soft object with pjpo 
= 0.1 [Figs. 2(a) and (d)] and the rigid object with p s /p 0 = 10, 
[Figs. 2(b) and 2(e)], the scattering amplitude \f(0)\ is circularly 
symmetric for k 0 a s = 0.1. This is due to the fact that scattering is 
dominated by the monopole term A 0 . The scattering behavior of soft 
and rigid objects is quite different for different values of k 0 a s . For the 
soft object [Figs. 2(a), 2(b), solid red lines] A 0 is still dominant for 
k 0 a s = 1, and no significant angular dependence is observed. For the 
rigid object, [Figs. 2(e), solid red line], higher order multipoles start 
to contribute quite significantly to the overall scattering when the size 
is increased to k 0 a s = 1. This has consequences on cloaking behavior. 
Consider the soft and rigid objects cloaked with shells with p c = 13 p 0 
and p c = — 2p 0 , respectively. For both shells, y = 0.5, i.e. a c = 2a s . 
Note that these values satisfy Eq. (10). Scattering amplitude \f(0) \ and 
normalized backward scattering amplitude \f(n)\/y/ch of the 
cloaked objects are shown in Fig. 2 (dashed blue lines). Scattering 
reduction is clearly observed for both objects in the long wavelength 
limit while in the Mie scattering regime only the soft object gets 
cloaked in every direction. 

Also, curves of \f (n) \j \fa s vs. k 0 a s in Figs. 2(c) and 2(f) show that 
the scattering from the rigid object undergoes a more dynamic vari- 
ation in /c 0 a 5 -range [0, 1], contrary to the scattering from the soft 
object, where only one broad maximum is observed in the same k 0 a s - 
range. This is associated with the fact that higher densities corre- 
spond to smaller effective wavelengths in the object, thus causing 
several internal resonances even for moderately sized objects. In 
addition, scattering reduction is more pronounced and broader for 
cloaking of the soft object [Fig. 2(c) and 2(f)]. 

We now analyze how introducing a suitably designed negative 
density shell can drastically lower the overall scattering. We schem- 
atically indicate the scattering, in analogy with the polarization in the 
electromagnetic case, with vertical arrows of opposite direction: the 
positive (solid) arrow stands for positive' scattered wave, whereas the 
negative (dashed) one corresponds to scattering with opposite phase 
[inset of Fig. 3(b)]. Figure 3(a) shows the contours of \pJPo\ (in 
logarithmic scale) that fulfill Eq. (10) for varying values of p s /p 0 
and y. The solid blue line in the figure represents the case p c = 0 
and + or - sign indicates whether the required shell density is positive 
or negative. Let e sca represent the scattering efficiency computed by 
the ratio of scattering cross-sections of the cloaked and bare objects. 
Figure 3(b) plots e sca vs. p c /p 0 for a cloaked object with p s /p 0 = 10 for 
various values of y and for k 0 a s = 0.1 (the quasistatic limit). Figure 3 
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Figure 2 | Scattering amplitude | in logarithmic scale for the soft object (pjp$ = 0.1) with (a) Aoa s = 0.1 and (b) Aoa s = 1 and for the rigid object 
(pJPo = 10) with (d) koa s = 0.1 and (e) koa s = 1. The dashed blue line represents the cloaked scenario whereas the solid red line stands for the bare 
object. Plots in (c) and (f) are the normalized backward scattering amplitude \f(n)\/y/a~ s vs. k 0 a s for soft and rigid objects, respectively. 



clearly shows that for certain values of p c , significant scattering 
reduction is achieved, and this is also obtained for different values 
of y. We notice that, for an ultrathin cloak (in the limit of y ~ 1), there 
is no dip in e sca ; instead the scattering reduction can be quite broad- 
band for thicker shells. For instance, the shell with y = 0.5 and p c /p 0 
= —2 significantly suppresses the scattering. This agrees well with 
the result given by Eq. (8), which considers only the zeroth-order 
multipole. As schematized in the inset of Fig. 3(a), a platonic cover 
with relative density less than unity may induce an out-of-phase 
monopolar mode with respect to the local elastic displacement field, 
thus permitting dramatic cancellation of the field scattered from the 



object. For other values of y, different values of p c /p 0 ranging from 
negative values to low positive (less than unity), are required to 
achieve significant cancellation of scattering. 

In reality, while modeling the frequency response of the cloaking 
phenomenon, one needs to take into account the inherent dispersion 
characteristics of the shell material since it is not possible to have a 
negative and frequency independent density due to causality con- 
siderations 47 . A negative density metamaterial may be achieved in 
acoustics by mimicking the way we achieve negative permittivity in 
electromagnetics with low frequency homogenization in high con- 
trast periodic media 47 " 52 , or high frequency 53 and high order homo- 




Figure 3 | (a) Contour plot of the solutions of Eq. (10). (b) Scattering efficiency e 3 " of the cloaked object with pjp 0 = 10 vs. pjp 0 for various values of y 
and k 0 a s = 0.1. The inset gives the sketch of cloaked object with arrows marking opposite phase of scattered wave in core and shell. 
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Wavenumber k Q a s 

Figure 4 | (a) The Drude-like density dispersion of the shell material with 
with pjpo = 10 vs. k 0 a s for various values of y e and y = 0.5. 

genization 54 in moderate contrast media. For elasticity, spring mass 
or inertial resonator metamaterials could lead to the desired effect 
near the resonance frequencies 55 . We suppose in the following that 
the material comprising the shell follows a classic Drude-like model 

K 

of the type p c (k 0 ) /p 0 = p in f — nnl — r [notice the fourth order 

dependence on k 0 stemming from the nature of the biharmonic 
equation], where p inf = 1, and k p and y e are the plasma wavenumber 
and Ohmic loss, respectively. This model can be obtained via a homo- 
genization applied to a composite shell with high contrast inclusions, 
and is in all ways analogous to that derived in 56 . This Drude disper- 
sion is plotted in Fig. 4(a) with parameters k p a s = 0.15 and 
y e = 10~ 2 kp. Figure 4 (b) plots scattering efficiency e sca of an object 
cloaked with a shell made of Drude material vs. k 0 a s for various 
values of y e . For this example, p s /p 0 = 10 and y = 0.5. Here, 
y e = 10~ 4 kp, (blue solid line), y e = 10 _2 /y (green dashed line), and 
y e = 10 -1 /y (orange dotted line) correspond to small, moderate, and 
high Ohmic losses, respectively. 

These curves are obtained using Eqs. (2) and (3), consistent with 
rigorous elastic scattering theory 46 and validated against results 



30 











Cloak ( Y =10" V) " 




— Cloak ( Y =10" 2 k 2 ) _ 




..-■Cloak ( y =10" 1 k 2 ) 




Cover 







0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

Wavenumber k a 
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15 and y e = 10~ 2 /y (b) Scattering efficiency e* ca of the cloaked object 



obtained using a finite elements commercial software 57 . It is clear 
that the presence of a cloak with negative effective density may allow 
for a drastic reduction of the overall scattering at a desired wave- 
number (k 0 a s = 0.11 in this example), independently of the angle of 
observation. It is stressed here that the scattering cross-section can be 
reduced by over five orders of magnitude compared to the uncloaked 
scenario and by six orders of magnitude compared to the case of an 
obstacle of the same size and density p c /p 0 = 10. A slight broadening 
of the cloaking dip and a corresponding deterioration of the effect 
may be noticed when elastic loss y e is increased. 

Discussion 

We have put forward a technique to cloak objects from elastic waves 
governed by the scalar fourth order biharmonic equation. The func- 
tionality of the obtained cloak could be seen in Fig. 5, where we plot 
the amplitude distribution of the scattered elastic displacement field 
U in the presence of cloaked [Fig. 5(a)] and uncloaked [Fig. 5(b)] 
objects. When it is surrounded by the negative density platonic shell 
(following the Drude-like model described above), the field ampli- 
tude is constant and equal to unity everywhere in space, in contrast to 
the considerably perturbed fields in the uncloaked case. The scatter- 



ed 



3= 

CD 



(b) 



k p a s = 0. 





(a) X (b) X 

Figure 5 | Time-averaged displacement field distributions (vector direction with white arrows, amplitude indicated by the contours) in linear scale of a 
rigid elastic object with pjpo =10 cloaked by a Drude-like shell (a) and on its own for comparison (b). The structures are illuminated with a unit- 
amplitude plane wave {k 0 a s = 0.11) propagating in the x-direction from right to left. 



SCIENTIFIC REPORTS | 4:4644 | DOI: 1 0.1 038/srep04644 



5 



ing reduction can be attributed to the proper choice of mass density 
of the platonic shell, in view of the scattering reduction predicted in 
Figs. 3 and 4. 

In summary, we proposed a design of an elastic cloak based on the 
scattering cancellation technique, inspired by earlier applications in 
electrodynamics and acoustic scenarios. Here, however, the chal- 
lenge is associated with the fact that a different, more complex ana- 
lytical form governs bending waves. The cloaking mechanism 
introduced here presents significant advantages in comparison with 
transformation acoustic designs: there is no need of anisotropy and 
inhomogeneity of the material parameters. Using a homogeneous 
isotropic platonic shell with low or negative density we have been 
able to greatly reduce the scattering from soft and highly rigid objects 
as well as from stress-free holes. It should also be noted here that the 
proposed technique cannot be used to cloak clamped objects in the 
quasistatic limit. However, cloaks made of layered shells have the 
potential to overcome this problem at higher frequencies (where 
scattering from clamped objects is finite) by adding extra degrees 
of freedom to cancel more than one scattering coefficient, as prev- 
iously demonstrated in electromagnetic cloak designs 58 . 

Experimental realization of this idea may be within reach in the 
near future (note for instance that lensing of bending waves via 
negative refraction was theoretically predicted using the biharmonic 
plate model in 41 ' 42 and experimentally confirmed in a thin 
Duralumium plate in 59 and a theoretical proposal for cloaking via 
geometric transforms in the biharmonic operator was also experi- 
mentally confirmed in 37,38 ), allowing for exciting applications in 
scenarios in which it is desirable to suppress the scattering from 
obstacles in thin plates in the motor vehicle and airplane industry, 
in which the scattering of bending waves in thin plates generated by 
engines may be totally suppressed, or smart secure buildings, pro- 
tecting them from bending and other types of seismic vibration 
damages caused by scattering from neighboring buildings 60 . 

Methods 

Analytical methods based on scattering Mie theory of cylindrical objects in thin elastic 
plates are used to obtain the numerical simulations in this Report (for further details 
and physical insight, see supplementary material 45 ). The vertical displacement of the 
plate is the solution of the fourth order partial differential equation of Kirchhoff. We 
proceed, as is usually done, by expanding the impinging plane waves and the scattered 
fields in terms of Bessel and Hankel functions in polar coordinate system centered 
with the object to be cloaked. We then apply four elastodynamic boundary conditions 
on each cylindrical interface in order to obtain the scattering coefficients for waves, 
which uniquely determine the displacement fields everywhere. The displacement 
field distributions and scattering cross sections are computed using Bessel develop- 
ments and Eq. (6) respectively. In the quasistatic limit, where the size of the elastic 
core sphere is much smaller than the wavelength and only the lowest- order Mie 
coefficient is important, an analytical formula is obtained [Eq. (10)] and it gives 
similar results to the full-wave simulations. Proper convergence for all the results is 
reached. The simulations given in Fig. 5 are obtained using the commercial finite- 
elements software COMSOL Multiphasics. 
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